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FOREWORD 


In Part 1 of this report (Schwiderski, 1978a), a combined hydrodynamical-empirical 
method was introduced to compute numerically harmonic partial tides in the world oceans 
with an accuracy of better than 5 cm, which is needed in various military and civil 
applications of today. In this report, the computed fortnightly lunar tide (Mf) is displayed 
in an atlas of tabulated tidal charts and plotted corange and cotidal maps. 
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ABSTRACT 


In Part I (Schwiderski, 1978a) of this report, a unique hydrodynamical interpolation 
technique was introduced, extensively tested, and evaluated in order to compute partial 
global ocean tides in great detail and with a high degree of accuracy. This novel method 
has been applied to construct the fortnightly lunar (Mf) ocean tide with about the same 
relative accuracy of better than 5 cm anywhere in the open oceans. The resulting tidal 
amplitudes and phases are tabulated on a 1° x 1° grid system in an atlas of 42° x 71° 
overlapping charts covering the whole oceanic globe. A corresponding atlas of global 
corange and cotidal maps is included to provide the reader with a quick general overview 
of the major tidal phenomena. The specifying hydrodynamical parameters of the model are 
listed along with quoted sources of empirical tide data, and significant tidal features 
are explained and discussed. As expected, the long-period fortnightly lunar tide (Mf) 
differs in structure considerably from the earlier computed semidiurnal and diumal ocean 
tides. In most areas, the oceanic tidal response resembles its forcing equilibrium tide, but 
amphidromic points do exist in the middle northen and southern latitudes. 


1. INTRODUCTION 


Part I of this report (Schwiderski, 1978a) introduced a unique combination of hydro- 
dynamical and empirical methods to model detailed ocean tides with a relative component 
accuracy of better than 5 cm anywhere in the open oceans. This enormous accuracy is 
well above minimum requirements set by, for instance, the National Aeronautics and Space 
Administration (NASA) and the Department of Defense (DoD) — to map the geoid at sea 
by satellite altimetry to within 10 cm. The following features of this unique hydrody- 
namical interpolation model made the achievement of this accuracy possible. 


a. A spherically graded 1° x 1° grid system is set up in connection with a cor- 
responding 1° x 1° bathymetry to assure a sufficient resolution of all important tidal 
phenomena. 


b. The bathymetry of the grid wise, simply connected ocean basin is hydrodynamically 
defined (Schwiderski, 1978c) by appropriate modifications of earlier realistic depth data ‘ 
collections. The hydrodynamical redefinition was needed in order to model the well-known 
strong distortion and retardation effects of shallow continental shelves, narrow ocean ridges { 
or island chains, and other significant bottom irregularities. ; 


c. The Boussinesq substitution of the turbulent Reynolds stresses is applied in the 
form of eddy dissipation with a novel physically meaningful eddy viscosity that depends 
linearly on the lateral grid-cell area and, hence, directly on the ocean depth. 


d. The linear law of bottom friction is introduced with a bottom-friction coefficient 
depending linearly on the bottom grid-cell area which is independent of the ocean depth. 
In boundary cells, the otherwise constant friction coefficient is subjected to an indirect 
cellwise adjustment in order to permit a consistent hydrodynamical interpolation (see h., 
below) of empirical tide data known from tide gauge stations at continental shores, 
islands, or other shallow-ocean bottom irregularities. 


e. The effects of the terrestial tide and the oceanic tidal load are included as simple 
second-order approximations in the sense of Love and Accad and Pekeris (1978). 


f. The Hansen-Zahel (Zahel, 1970 and 1977; Estes, 1977) finite differencing technique 
is modified by a new differencing scheme in time which improves decay, dispersion, 
and stability characteristics of the numerical procedure and facilitates the simple indirect 
adjustment of the bottonrfriction coefficient in the hydrodynamical interpolation technique 
(see d. and h.). 
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g. At land-ocean cell walls, the conditions of no-flow across and free-slip along the 
boundaries are enforced. The no-flow condition is subsequently relaxed by allowing con- 
trolled periodic inflows and outflows over the mathematically assumed boundaries. This 
allowance redefines indirectly more realistic shorelines in order to further improve the 
consistency of the hydrodynamical interpolation of empirical data (see d. and h.). 


h. A unique hydrodynamical interpolation technique is introduced which incorporates 
into the theoretical model empirical tidal constants collected from over 2 000 tide-gauge 
stations around the world in a hydrodynamically consistent fashion (see d., f., and g., above). 


i. A new higher-order approximation of Arctic Ocean tides is used, that is described 
in Schwiderski (1981h). 


With these features, the new model was successfully applied to chart the semidiurnal 
principal lunar (Mz) ocean tide with the desired accuracy. The technique and accuracy of 
the model were extensively described and discussed in Part I of this report as well as in 
subsequent publications and symposia presentations by the author (Schwiderski, 1978a, b; 
1979a, b, c, d, e; and 1980). 


Essentially the same hydrodynamical interpolation technique has been applied to chart 
the fortnightly lunar (Mf) ocean tide (see Section 2). However, it must be emphasized 
that the accuracy of the Mf model is less verified by empirical tide data than the accuracy 
of the earlier computed semidiurnal and diurnal partial tides (Schwiderski, 1979e; 1981a-g). 
In fact, reliable empirical tide data were only available for relatively few tide gauge 
stations. Nevertheless, particularly in the open Pacific, the accuracy of the Mf model is 
very well substantiated by the recent empirical tide data published by Luther (1980) for 
widely scattered island stations. The thorough harmonic analysis of Luther improved older 
tide data in some cases significantly and, hence, was most valuable in the computation of 
the tide model. 


A complete listing of all suerces of empirical ocean tide data, which were interpolated 
into the Mf tidal charts, is presented in Appendix A. In the meantime, Section 2 of 
this report lists the significant hydrodynamical input parameters that specified the con- 
structed Mf ocean tide. The major features of the global Mf tide are discussed in Section 3. 
A complete numerical display is presented in Appendix A where all tidal amplitudes and 
phases are gridwise tabulated in map-like charts. Corange (equi-amplitude) and cotidal 
(equi-phase) maps of the Mf ocean tide are plotted in Appendix B. 


2. Mf OCEAN-TIDE PARAMETERS 


The astronomical fortnightly lunar (Mf) equilibrium tide 1 (or tide-generating potential 
Gn; see Schwiderski, 1978a) at the geographical point (A, ¢) and instant (Y, D, t) is 
determined by 


n= K (5 eU cos (of + x) (1) 


= 9.81 m/sec? earth gravity acceleration 
= longitude (east in rad) 
= latitude (north in rad) 
(= 1975) = year number 
= day number of year Y(D = 1 for January 1) 
= universal standard time of day D (in sec) 
= 0.041 742m = Mf equilibrium tide amplitude 
0.53234°10-5 sec"! = Mf tide frequency 
2180 / 180 = Mf astronomical argument (in rad) 
270.434 358 + 481 267.883 141 377 - 0.001 13372 + 1.9 . 106 73 
mean longitude of the moon relative to Greenwich midnight of day D (in deg) 
T = [27 392.500 528 + 1.000 000 035 6D] 1/36 525 
D = D 365 (V 1975) + Int IT- 1973)/4] 
Int[x] = integral part of x 
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The corresponding instantaneous ocean partial tide (Schwiderski, 1978a) is determined 
by 


g = & cos (ot + x - 5) (2) 
where the local harmonic constants 
t = (A, ) = Mf ocean tide amplitude (in m) 


and 


an 


= 6(A, G) = Mf ocean tide Greenwich phase (in rad) 


must be determined, say, by linear interpolation in the tidal charts of Appendix A. 


A simple second-order approximation in the sense of Love and Accad and Pekeris 
(see Part I, Schwiderski, 1978a, 1979c, and 1980; and Accad and Pekeris, 1978) yields 


& =~ 0.612n and e = -0.0667§ (3) 


i. e., the corresponding terrestrial tide (© and the earth dip (°° (yielding) under the oceanic 
tidal load §, respectively. A more elaborate and probably slightly more accurate earth dip 
may be computed by using Farrell’s Green function (see Farrell, 1972 and 1973; 
and Schwiderski, 1980). In linear superposition, one finds the corresponding instantaneous 
geocentric partial Mf tide: 


e (4) 


A detailed description of the hydrodynamicalempirical model to compute the ocean 
tide amplitudes ? and phases 6 (listed in Appendix A) was given in Schwiderski (1978a, 
1979c, d, and 1980). In particular, all model input parameters such as the dimensionless 
eddy coefficient e (Eq’s. 103 and 123), the bottom-friction parameter b (Eq’s. 4a and b), 
and the differencing parameters x and & (Eq’s. 64 and 72) were all specified in Schwiderski 
(1978a) (referenced equations). These parameters were determined for M, by extensive 
trial-and-error computations and zemained unchanged for the construction of all semidiurnal 
(M,, Sz, Nz, Kz) and diurnal (K,, O1 Pi, Qi) tides (see Schwiderski, 1978a, 1979e, 
1981a-g). However, for the computation of the much slower long-period tide Mf, the 
eddy and bottom friction coefficients had to be reduced by 33% and 50%, respectively. 


In the computation of the Mf tide model, the following mode-dependent parameters 
were used (see referenced equations in Schwiderski, 1978a): 


a. The time step At (Eq’s. 64, 123) 
At = 182.1444 sec (5) 


b. The hydrodynamical interpolation control limits, k,, k,, and k; (Ed's. 88, 89, 
94, 97, and 99) 


k, = 0.010, k, = 0.020, & = 0.5 (6) 


It may be noted that the best“ fit of the model to the empirical tide data required a 
considerable reduction of the eddy and bottom friction coefficients and of the most 
important control parameters k, and k,. While all semidiurnal tide models required 
k, = 0.045 and k, = 0.045, the diurnal tides needed a first reduction to k, = 0.025 and 


k, = 0.040, with the additional reduction shown by Ed. (6) for Mf. As was mentioned 
for the K, tide (see Schwiderski, 1981b), a reduction of this sort was expected since 
these parameters characterize the intensity of the turbulence (eddy dissipation) which is 
present in all tidal motions (see Schwiderski, 1978a, 1979c, d, and 1980). As is physically 
well known, the activity of turbulence depends on the velocity of a fluid motion. When 
the fluid velocity is significantly accelerated, as in tidal currents from long-period to 
diurnal and to semidiurnal tides, the vigor of turbulence increases accordingly. 


During the construction of the Mf tide model, a major computational difficulty 
became apparent. The accurate determination of the amplitudes and phases of all semi- 
diurnal and diurnal partial tides required (see Schwiderski, 1978a, 1979c, and 1980) two 
accurate time-dependent ocean tidal states about one-quarter period apart. To compute 
such final ocean states for the fortnighly Mf tide with a time step of At ~ 3 min (Eq. 5) 
would have required an unacceptable computer time. The planned computation of the 
; monthly moon (Mm) tide and semiannual sun (Ssa) tide would certainly be infeasible. 
„ To overcome this difficulty, it was found far less costly to compute two suitable ocean 
tidal states independently from two initial states of one-quarter period difference. Ad- : 
ditional averaging procedures and tidal flipping by one-half period (multiplication by 
1) were advantageously applied to speed up the convergence to the necessary final „ 
periodic state. 
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3. Mf OCEAN-TIDE FEATURES 


The entire constructed Mf ocean tide is gridwise displayed in map-like amplitude and 
phase tables in Appendix A. The 42° x 71° charts cover the whole globe north of 
colatitude 169° (Antarctica) in three zones: a northern zone N from 0° to 71° colatitude, 
a middle zone M from 48° to 118° colatitude, and a southern zone S from 98° to 168° 
colatitude. The overlapping geographical areas of the tidal charts have been chosen to 
provide a worldwide coverage for special applications and to allow the reader to scan the 
large amplitude and phase charts together in order to evaluate their quality and visualize 
the important tidal features. In addition, a generally superficial overview of some tidal 
features can be recognized by inspecting the more schematically plotted corange and 
cotidal maps provided in Appendix B. 


The tidal charts and maps permit the viewer to follow the tidal waves, that is, the 
high water fronts (crests), in forward (or backward) direction, for instance, on their 
rotation around the amphidromic points. In the tidal phase charts of Appendix A, it is 
best to start from the prominently visible 0° = 360° or 100° cotidal lines. Since the 
Greenwich phases specify the time lags (in degrees: 30° ~ 27 hours) of the tidal crests 
relative to the cresting time of the corresponding equilibrium tide over the central latitudes, 
one gathers a vivid impression of the significant global and local tidal phenomena. 


For an easy evaluation of the tidal charts in Appendix A, all hydrodynamically 
interpolated em irical tidal amplitudes and phases have been visibly marked by subbars. 
The approximate geographical locations of the important amphidromic points of zero 
amplitudes are marked by a circled (x). By following the tidal waves on their periodic 
rotations, one finds these waves passing through the specially marked stations in empirically 
correct time and with the correct height. In fact, all empirical tidal phases and amplitudes 
are coherently integrated. This is particularly impressive for the charts of the Pacific Ocean, 
where the improved empirical data of Luther (1980) from many scattered island stations 
fit smoothly into the surrounding computed tides. From the smoothness features of er- 
ratically interpolated tidal data (see Parts I and II), one concludes that this result is 
not an artifact of the interpolation applied but constitutes a vivid manifestation of the 
excellent compatibility of both the empirical and hydrodynamical procedures combined. 


The corange and cotidal maps of Appendix B may be compared with the corre- 
sponding maps of Kagan et al. (1976) and Agnew and Farrell (1978). While there is some 
qualitative agreement with both models, the present maps clearly confirm the conclusions 
of Luther (1980). Particularly, over the central oceans, the Mf ocean tide differs con- 
siderably from its driving equilibrium tide (Eq. 1), which is simply a standing wave with 


the nodal lines (amphidromic lines) along the latitudes G ~ 35.26°. In the neighborhood 
of these nodal lines, the Mf ocean tide is very low (1 mm or less) and the phases 
appeared, of course, much less determined than in higher amplitude areas. Nevertheless, 
the existence of several amphidromic points near the nodal lines, particularly in the 
southern oceans, ren: ained quite steady during the final steps of the computations. 
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APPENDIX A 


ATLAS OF 1° x 1° Mf OCEAN-TIDE AMPLITUDE 
ö AND PHASE CHARTS FOR 42° x 71° AREAS 


1. GUIDE TO TIDAL CHARTS 


=  m:Longitude Number 

=  n:Colatitude Number 

= (m- 0.5): Geographical Longitude East 

= (n - 0.5)°: Geographical Colatitude 

EA,,> 6): Amplitude (in mm) 

=  6(A,,, 9): Greenwich Phases (in deg; 30° 27.3h) 


2 


Amphidromic Points 
=  Subbars Mark Empirical Input Data at Shore Stations 


1 ® gt Ze 
2 
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2. SOURCES OF EMPIRICAL TIDE DATA 


Wunsch (1967), International Hydrographic Bureau (1978), and Luther (1980). 
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Amplitudes k of corange lines in mm. 


Greenwich phase & of cotidal lines in 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 
165, 180, 195, 210, 225, 240, 255, 270, 285, 300, 315, 330, 345, 360 = 0 where 
30° 27.3 hours. 
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